Abstract. This is the first of two papers in which we ad- 
dress the physics of the low-frequency Quasi-Periodic Os- 
cillation (QPO) of X-ray binaries, in particular those host- 
ing a black hole. We discuss and repeat the recent analysis 
and spectral modelling of the micro-quasar GRO J1655- 
40 by Sobczak et at (2000, hereafter SMR), and compare 
it with GRS 1915+105; this leads us to confirm and ana- 
lyze in more detail the different behavior noted by SMR, 
between GRO J1655-40 and other sources, when compar- 
^vi ing the correlation between the QPO frequency and the 
. disk inner radius. In a companion paper (Varniere et al. 
' , 2002, hereafter Paper II) we will show that these op- 
s ! , posite behaviors can be explained in the context of the 
?H ' Accretion-Ejection Instability recently presented by Tag- 
ger and Pellat (1999). 

We thus propose that the difference between GRO J 1655- 
40 and other sources comes from the fact that in the for- 
£S| ' mer, observed in a very high state, the disk inner radius 
always stays close to the Last Stable Orbit. 
In the course of this analysis, we also indicate interesting 
, differences between the source properties, when the spec- 
^v^j ' tral fits give an anomalously low inner disk radius. This 
£C} , might indicate the presence of a spiral shock or a hot point 
C*~> ' in the disk. 

o . 

£^ ' Key words: Accretion discs - X-rays : binaries - Stars 
^~ ! individual GRS 1915+105; GRO J1655-40 
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1. Introduction 

The observation of Quasi-Periodic Oscillations (QPOs) in 
X-ray binaries is widely considered as a key to under- 
standing the physics of the inner region of accretion disks 
around compact objects. Although the origin of the QPOs 
is still debated, it is considered that such models should 
yield information of prime importance on the physics and 
the geometry of the accretion process, in these and other 
accreting sources. In particular, in the last years, many 
observational results have pointed to the importance of 
a low-frequency QPO, found in binaries hosting neutron 
stars or black holes. A recent result, reported by Psaltis et 
al. (1999), shows that in these sources its frequency is cor- 
related with that of a higher-frequency QPO (in neutron- 
star binaries, this is the lower of the pair of so-called 'kHz 
QPOs'), believed to be close to the rotation frequency at 
the inner edge of the disk. Thus the frequency of this QPO 
ranges from ^ 1 to a few tens of Hz in black-hole and 
neutron-star binaries. It has been intensively studied, in 
particular, in the micro-quasar GRS 1915+105, where its 
occurrence in the low-hard state is so frequent that it has 
been dubbed the "ubiquitous" QPO (Swank et al. , 1997; 
Markwardt et al. , 1999, hereafter SM97 for both these ref- 
erences). These studies have allowed SM97 and Muno et 
al. (1999) to find a number of correlations, showing that 
the QPO probably has its origin in the disk, although it 
affects more strongly the coronal emission. 

One of these correlations involves the QPO frequency 
and the color radius of the disk r co i , considered as a mea- 
sure of the disk inner radius Ti nt . Even though this is not 
a very tight correlation, it would mean that the QPO fre- 
quency would most often correspond to a Kcplerian ro- 
tation frequency in the disk, at a radius of a few times 
(typically ~ 2 to ~ 5) ri n t] this is consistent with the cor- 
relation of Psaltis et al. (1999), where the high-frequency 
to low-frequency ratio is about 11, assuming, as usually 
believed, that the high-frequency QPO is close to the Ke- 
plerian frequency at the inner edge of the disk. 
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However these observational results must be taken 
with extreme caution: r co i is obtained from a fit of the 
source spectrum with a model, and the results of the fit 
raise many questions: 

1. The model spectrum consists of a multi-color black- 
body part, thought to originate in an optically thick 
disk, and a power-law tail at higher energies, consid- 
ered as inverse Compton emission from a hot corona 
which might lie either above the disk, or in the re- 
gion between the disk and the central object. We will 
call them, for simplicity, the disk and coronal emis- 
sions respectively. However in the low-hard state where 
the low-frequency QPO is most often seen, the coro- 
nal emission often dominates the disk one, even at low 
energies. In this case the extraction of the disk param- 
eters becomes difficult. 

2. Indeed the fits often give values of r co i which must 
clearly be ruled out, since they are extremely low and 
sometimes even within the Schwarzschild radius. 

3. In a recent paper Merloni et al. (2000) (hereafter MFR) 
have shown that other effects, associated with the ver- 
tical structure of the disk, could strongly affect the 
determination of the disk parameters. Using synthetic 
spectra to test the fitting algorithm, they find that the 
disk black body model underestimates the radius by 
a factor which may be quite large, when the coronal 
emission dominates the disk emission. This makes the 
extraction of disk parameters difficult, since we have 
little information on the conditions causing this effect. 

4. Even when the disk emission is strong enough to allow 
a reliable fit, the model spectrum contains assump- 
tions which may affect the outcome. In particular the 
disk is assumed to be axisymmetric with a tempera- 
ture varying as T -3 / 4 ^ following the standard ev-disk 
model. These are strong assumptions which may af- 
fect the validity of the fit parameters. However, since 
the RXTE/PCA data starts above 2 kcV, the obser- 
vations are not very sensitive to the whole structure of 
the disk itself. On the other hand, various corrections 
(e.g. due to electron scattering) must be applied. One 
can only expect that, when the assumptions are not 
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too far from reality, there is a proportionality between 
the color radius r co i and the disk inner radius ri nt . 

Points (1-3) above thus lead us to take with extreme 
caution, or rule out, determinations of the disk parame- 
ters when the ratio of the disk black body flux (determined 
from the fit) is less than about half of the total flux. We 
will see below that indeed, among the data points we ob- 
tain for GRO J 1655-40, all the ones giving anomalously 
low values of r co ; are excluded if this criterion is applied. 
On the other hand point (4) above is much more difficult 
to address, as long as we have no observational tests of the 
disk structure. Determining the nature of the QPO would 
be very likely to shed light on this, but this also depends 
on our knowledge of the disk parameters. 

Recently SMR have analyzed observations of the 
micro-quasar GRO J1655-40, and found that in this source 
the low-frequency QPO differed markedly from other 
sources, in particular XTE J1550-564: the color radius was 
found to be quite small (down to a few kms), and to in- 
crease with the QPO frequency. On the contrary in XTE 
J1550-564, and in most observations of other sources, the 
color radius and the QPO frequency are inversely corre- 
lated. 

In this paper our goal is to check this reversed corre- 
lation which (as will be seen in Paper II) might be consis- 
tent with the explanation of this QPO by the Accretion- 
Ejection Instability (AEI) recently presented in Tagger & 
Pcllat (1999). It would be due to relativistic effects, when 
the disk inner radius approaches the Last Stable Orbit. 



erences therein). We have thus analyzed two such cycles, 
essentially repeating the analysis of SM97 and turning our 
attention to the frequency-radius correlation. 

We will present the observations of GRO J 1655-40 and 
GRS 1915+105 in sections O and O respectively. We 



Thus we begin in section 2.1 by reconsidering the data 



points from GRO J1655-40 presented by SMR, repeating 
the whole data reduction and fitting process; we discuss 
them in regard to the criterion of MFR, which we take as 
a guide for the validity of the fitting procedure. This re- 
sults in excluding all the points at anomalously low radius 
values, and retaining points only in a limited radial range. 
We then turn to the other source studied by SMR, XTE 
J 1550-564, but the criterion of MFR leads us to exclude 
almost the whole data set: the source was at the time of 
the study in an extremely bright high soft state/ very high 
outburst state. We thus chose to study a different source, 
GRS 1915+105, because its high variability and the large 
number of observations including the low frequency QPO 
suggested that we could explore the radius-frequency cor- 
relation over a much broader range of values for r co ;. In 
particular the observations of SM97 showed that, during 
the much studied 30 min. cycle of that source, both vqpo 
(the QPO frequency) and r co i vary by a large factor, and 
apparently in opposite directions. The choice of this par- 
ticular observation was driven also by the existence of 
a multiwavelength campaign (Mirabel et al., 1998), that 
could allow a richer discussion between QPOs, disk insta- 
bilities and ejections of plasma, since several authors have 
now linked the changes of luminosity to disk instabilities 
(e.g. Fcroci et al. , 1999; Fender et al. , 2001, and ref- 



will then in section 2.3 discuss and compare the results, 



which will be then be confronted with a theoretical model 
in Paper II. We will also discuss the different behaviors 
observed when the disk color radius takes anomalously 
small values, leading to a discussion of disk properties and 
to a possible explanation for these very small values. 



2. Data reduction and analysis 

2.1. GRO J1655-40 

GRO J1655-40 (Nova Scorpii 1994) was discovered in 1994 
by the Burst and Transient Source Experiment (BATSE) 
(Wilson et al. 1994). Its optical counterpart was discov- 
ered in 1995 (Bailyn et al.,1995), and the orbital param- 
eters gave a primary mass of 7.02 ± 0.22M Q (Orosz & 
Bailyn, 1997; van der Hooft et al., 1998), placing it well 
above the 3M Q maximum mass for neutron stars. Its dis- 
tance is estimated at 3.2 ± 0.2 kpc (Hjellming & Rupen, 
1995), with an inclination angle - 69.5°. GRO J1655-40 is 
also known as one of the few galactic microquasars, since 
radio observations show two extended radio lobes with su- 
pcrluminal motions (Hjellming & Rupen, 1995). In order 
to be able to compare directly the results on the disk color 
radius with the ones presented below for GRS 1915+105, 
we have re-done the data reduction and analysis of the se- 
ries of twelve public observations (AOl) of GRO J1655-40 
during its 1996 outburst (see figure Q), already presented 
by SMR, when QPO were observed. Figure |l| shows the 
BATSE and RXTE ASM light curves during that period, 
and the dates where QPO were observed. Table [l] gives 
the corresponding disk parameters. 

For our analysis the RXTE spectra were obtained us- 
ing only the Proportional Counter Array (hereafter PCA). 
We used the standard 2 data only when the offset point- 
ing was less than 0.02° and the elevation angle above 10°; 
in almost all cases all the PCU were "on" during the ob- 
servation, but we chose to use only the spectra extracted 
from PCU and 1 separately, and then fit them together in 
XSPEC. Given the high flux of the source this allows us to 
get sufficient statistics while limiting the effects associated 
with the various PCUs. The data reduction/extraction 
was made following the standard RXTE data reduction. 
In order to compare precisely with the results of SMR, we 
have used FTOOLS package 4.1, but checked that package 
5.0 gave nearly identical results. 

We first generated spectra using all the PCU layers, 
which gives a higher count rate, and then extracted spec- 
tra only from the top layer; the fits parameters deduced 
with the two methods presented no significant differences, 
and since there was, for this source, no particular need 
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Fig. 2. Spectral fits for the first two observations of GRO J1655-40. On July 27 the black-body part dominates 
at low energy, allowing a satisfying extraction of disk parameters. On August 1 the coronal part dominates at all 
energies. The plots show (top curve) a superposition of spectra from PCU's & 1, and (lower curves) the spectral 
components (multicolor black body and power-law tail) obtained from the fit. The X axis is in keV, and the Y axis is 
in Counts/sec/keV. 




10200 10300 10100 10500 



Truncated Julian Days TJD 10083 = 01 01 1996 

Fig. 1. GRO J1655-40 light curves during 1996, obtained 
at high energies with BATSE (20 - 100 keV) and at lower 
ones with the RXTE All Sky Monitor (1.2-12 keV). 
Arrows show the dates where low-frequency QPO were 
observed. 

to get the highest count rate (the source was sufficiently 
bright, and the exposure sufficiently long that the data 
had small statistical errors), the results presented here are 
obtained with the later method. The background was esti- 
mated using the PCABACKEST tool version 2.1, and we 
chose to ignore channel energies bellow 2.5 keV , and above 
20 keV, in order to use the best detector energy band 
(choosing this band, we avoided the uncertainties due to 
the response matrix at lower and higher energies, and the 
problems due to background subtraction near 25 keV - 



see also SMR); we then added a systematic error of 1.5% 
on each data point. The response matrices were generated 
using the PCARSP tool version 2.38, and we deadtime cor- 
rect the spectra following the RXTE cook book method. 
The spectral fits were obtained using XSPEC version 10., 
using the multicolor disk blackbody model with a power 
law tail at higher energy, and a hydrogen column absorp- 
tion fixed at N H = 0.9 x 10 22 cm -2 as found with ASCA 
(Zhang ct al., 1997a). We have found that no other con- 
tributions to the spectrum (e.g iron line) needed to be in- 
troduced in the fit, giving results consistent with those of 
SMR while reducing the number of fit parameters. On the 
other hand, since there is little uncertainty on the QPO 
centroid frequencies, we have not re-done the power spec- 
tra analysis, and we use here the frequency values given 
by SMR. 

Table [l] shows the results of our fits. Column 7 shows 
the ratio (Disk black body flux/total flux) which, accord- 
ing to MFR, we will first use to discriminate between two 
different spectral states of the source. For our purpose, 
the total flux is obtained by extrapolating the spectra to 
50 keV. Table [l] shows that the observations clearly come 
in two categories: one where the ratio Disk black body 
flux/total flux is lower than .5, so that according to the 
criterion of MFR the disk radius value should be consid- 
ered as unreliable, and where indeed the radius is always 
unrealistically low; and one where the ratio is higher than 
.5 and where the obtained radius appears to be more trust- 
worthy: it lies in the range we expect for a disk around a 
black hole, and gives more consistent values. 

Figure ^ shows the spectra for the first two obser- 
vations, typical of these two behaviors: on July 25 the 
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disk emission dominates at low energy, allowing a sat- 
isfying spectral fit. On August 1 the coronal emission 
dominates at all energies, and results in an unreliable 
fit. Figure || is the resulting plot, for the data points we 
retain, of the QPO centroid frequency versus the disk 
color radius, where r co i is given in km, with R* = 
((r coi I km) /(D/Wkpc)) 2 cos 6 R* being the XSPEC di- 
mensionless parameter. Comparing with figure 2 of SMR, 
we note that the points we exclude correspond to higher 
temperatures and lower disk flux, ft is quite remarkable 
that they are also the observations where a higher fre- 
quency QPO is simultaneously observed. This might con- 
firm that indeed the disk is in a different state in the obser- 
vations we retain and the ones we reject, but we will not 
attempt to explore an explanation, which would clearly 
require reaching a much better understanding of the disk 
behavior than our present one - although in our discussion 
we will contrast this with the behavior of GRS 1915+105. 
On the other hand, with the data points we retain, we 
confirm the general tendency noted by SMR that, con- 
trary to other sources, the QPO frequency increases with 
increasing disk radius in GRO J 1655-40. 

Fig. 3. Plot of vQPO vs r col for GRO J1655-40. X axis is 
in km, Y axis is in Hz. 



2.2. GRS 1915+105 

We thus turn to the microquasar GRS 1915+105. This 
source has been the object of very frequent observations, 
many of them showing the low-frequency QPO and strong 
variability. In particular, the very exhaustive analysis of 
Muno et al. (1999) shows the QPO, in different spectral 
states, and over a large range of frequencies. Furthermore, 
SM97 had shown that, during the ~ 30 minutes cycles of 
GRS 1915+105, the QPO could be tracked at a varying 
frequency, together with a color radius varying by a factor 
~ 5, during the "low and hard" part of the cycle. These 
cycles have been extensively studied; multiwavelength ob- 
servations (Mirabel et al., 1998; see also Chaty, 1998) have 
shown them to be associated with supra-luminal ejections, 
and they turn out to be the most frequent type of behav- 
ior of that source after the steady ones (Muno et al, 1999; 
Belloni et al, 2000). Observing the correlated changes of 
radius and QPO frequency during such events may limit 
the influence of other effects, such as the mass accretion 
rate etc., and provide a way to analyze the disk history 
during such events. This could then complement the anal- 
ysis of Belloni et al. (2000), based only on color-color di- 
agrams. 

Figure ^| shows three consecutive cycles observed by 
RXTE on 1997 Sept. 9. The first cycle is the one shown in 
detail by SM97, while the second cycle was the object of 
the multiwavelength observations of Mirabel et al. (1998). 



We have thus chosen to analyze these two consecutive cy- 
cles: we repeat the analysis of SM97 for the first cycle, 
using for the spectral fits and color radius derivation the 
same procedure as for GRO J1655-40. In order to con- 
firm the results we obtain, we repeat the same analysis 
for the following cycle. We have also made a rapid analy- 
sis showing that the same behavior is observed during the 
third cycle. We treat only the time interval over which 
the QPO was observed by SM97, i.e. starting near the 
transition from the high/soft to the low/hard state, and 
ending at the intermediate peak halfway through the low 
state. At that time the QPO disappears, and the power- 
law emission decreases dramatically (Mirabel et al., 1998; 
Chaty, 1998). 
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Fig. 4. Lightcurve of GRS 1915+105 during the whole 
observation of 1997 Sept. 9. X axis is in seconds since the 
beginning of the observation, Y range is in count/s. 



2.2.1. Spectral analysis 

Our goal was to track the variations of the QPO fre- 
quency and of the disk color radius, while they vary rather 
rapidly. We thus reduced the time interval over which the 
spectral fit and the QPO frequency determination were 
obtained, to the minimal value compatible with sufficient 
statistics. We used the standard 2 data accumulated over 
16s intervals; although the source occasionally varies on 
shorter time scales, we have checked that this does not 
affect the results we present. Due to the poorness of the 
statistics, we extracted spectra from all available PCU 
(though only PCU to 3 were "on" during the whole 
observation), and all layers to get the most possible 
incoming flux. The spectra, after dead time correction 
and background correction, were first fitted, as for GRO 
J1655-40, between 2.5 and 20 kcV in XSPEC with the 
standard model consisting of a multicolor disk blackbody 
and a power law tail at higher energy, taking into account 
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# 


Date 


T col (keV) 


Rcol 




Powerlaw 


Ratio 


Xi (91 d.o.f) 


1 


07-25-96 


1.272 ± 3.6 x 10~ 3 


20.62 ±0.15 


2.589 ±0.022 


15.93 ±0.98 


0.64 


0.69 


2 


08-01-96 


1.63 ± 1.4 x 10" 2 


8.41 ±0.28 


2.49 ±0.012 


50.68 ± 1.78 


0.21 


0.76 


3 


08-06-96 


1.53 ±8.7 x 10~ 3 


10.90 ±0.23 


2.65 ±0.013 


47.93 ± 1.856 


0.31 


0.7 


4 


08-15-96 


1.247 ± 3.7 x 10 -3 


20.69 ±0.15 


2.56 ±0.024 


13.58 ±0.9 


0.65 


0.98 


5 


08-16-96 


1.27 ±6.6 x 10~ 3 


16.88 ±0.27 


2.59 ±0.016 


34.40 ± 1.578 


0.36 


0.82 


6 


08-22-96 


1.25 ±5 x 10~ 3 


19.10 ±0.20 


2.6 ±0.017 


26.79 ± 1.26 


0.47 


0.62 


7 


08-29-96 


1.59 ± 1.2 x 10~ 2 


9.41 ±0.27 


2.54 ±0.013 


50.28 ± 1.86 


0.24 


0.68 


8 


09-04-96 


1.25 ±4 x 10" 3 


20.56 ±0.17 


2.61 ±0.02 


19.50 ± 1.0 


0.59 


0.72 


9 


09-20-96 


1.23 ±3.6 x 10" 3 


21.10 ±0.14 


2.49 ± 0.023 


9.94 ±0.65 


0.69 


0.7 


10 


10-15-96 


1.29 ±3.5 x 10~ 3 


20.80 ±0.13 


2.36 ±0.028 


7.098 ± 0.56 


0.74 


0.6 


11 


10-22-96 


1.28 ±3.5 x 10~ 3 


21.13 ±0.13 


2.48 ± 0.025 


10.79 ± 0.76 


0.70 


0.62 


12 


10-27-96 


1.32 ±5 x 10" 3 


17.84 ±0.21 


2.64 ±0.017 


34.31 ± 1.66 


0.46 


0.7 



Table 1. Best fit parameters for the observations of GRO J1655-40. R* = ((R co i / km) / (D /lOkpc)) 2 cos 6, where 9 is 
the inclination angle. Powerlaw normalization is in unit of photon/ cm 2 /s at 1 keV. 



a low energy absorption with Njj fixed at 5.7 x 10 22 cra -2 
( Markwardt et ai, 1999). The fits were not satisfactory, 
with an average reduced x square ~ 2.5. An analysis of 
the residual showed a large difference between the spectra 
and the model used around 6 keV, so that in a second 
pass we added to the model an iron emission line, with 
centroid energy found at ~ 5.9 keV. 

As for GRO J1655-40, the coronal emission dominates 
the light curve in some observations. We thus choose the 
following procedure: we distinguish in the observations a 
high state and a low state; the latter is defined by the 
combination of a low temperature ( 1 keV) and a large 
radius. In this state, the coronal contribution dominates 
but the correction of MFR (which was derived for high 
states) should not apply, and the color radius is much 
larger than at other times, so that we retain these points - 
which correspond to most of the intervals where the QPO 
was present. In the high state on the other hand, the ra- 
dius is much lower and we apply the same criterion as 
above, retaining only data points where the black-body 
contribution exceeds half of the total (extrapolated to 50 
keV), although we find that in general the radius values 
for the rejected points are more consistent than in GRO 
J1655-40. The resulting parameters and classification are 
presented in Table |[ 

2.2.2. Temporal analysis 

As we wanted to follow the QPO in the range 2-12 Hz 
found in SM97, we generate power spectra for the first cy- 
cle of figure U using the binned data, with 8 ms resolution, 
and the event data with time resolution of 62 {is collected 
during 4s ccumulation, after background correction. Ac- 
tually the background is so low that retaining it would 
change the results only by less than ~ .5%. From this 
we generated a dynamical power spectrum of the source, 
repeating the analysis of SM97. Our result, shown on fig- 
ure H (left), confirms theirs: the shaded, U-shaped pattern 



seen on the lower part of figure || between times ~ 1000 
and ~ 1500, shows the QPO frequency and its evolution 
during the cycle. We then repeated this analysis for the 
second cycle, confirming the generality of this behavior 
(figure ||, right). 

The second step was then to track as precisely as pos- 
sible the QPO frequency, and to correlate it with the disk 
color radius extracted with the same process used for GRO 
J1655-40 in section 2.1. We do this using the same time 
resolution data, but with exposures of 16s. Power spec- 
tra were then generated using POWSPEC version 1, after 
background subtraction and collimator correction; they 
were then fitted with a model consisting of a power law 
with average index «~ —1.5, plus a lorentzian whose line 
center is taken as the QPO centroid frequency; in some 
cases several lorentzians were required to fit the 16s power 
spectrum; we then chose the one closest to the frequency 
at which power was strong in power spectra of neighboring 
intervals. The trend of the usually strongest feature is eas- 
ily picked out from the dynamical power spectrum. The 
choice of this 16 second interval for the temporal analysis 
is a compromise between conflicting constraints: looking 
at variations of the disk radius and the QPO centroid fre- 
quency on short enough time scales (without averaging 
on too long time intervals over which the source varies 
rapidly), having enough flux to get sufficient statistics on 
each data point, and getting the QPO frequency with good 
accuracy. 

2.3. Results 

The individual power spectra, obtained from interval 
# 2, are shown in Fig. [| One sees that, with the time 
intervals of 16s., the QPO can most often be identified 
and properly found with our fitting procedure. Figure 
[| shows superpositions of our results for the disk color 
radius, the QPO frequency, and the disk temperature 
during both cycles. Figure || shows, for both cycles, the 
correlation between the QPO frequency and the disk color 
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Fig. 5. Dynamical power spectra of GRS 1915+105. The left and right panels correspond respectively to the first 
and second cycles analyzed. The X axis is is seconds, with the zero corresponding in each case to the beginning of 
accumulation. 



radius. At large r co i the QPO frequency decreases with 
increasing r co i, as discussed by Sobczak et al. (2000) for 
XTE J1550-564; this is also the general tendency for this 
correlation, noted by SM97 during one cycle, and more 
generally by Muno et al. for an ensemble of observations 
(in particular for similar 30 min. cycles). 

It would be tempting to interpret the points at lower 
radius as a turnover, indicating a reversed correlation as 
seen in GRO J1655-40. However, returning to Fig. [t], one 
sees that this interpretation must be excluded: the points 
at lowest radius correspond to the very first manifesta- 
tions of the QPO. At this time the count rate has not yet 
started to decrease (this might indicate that the appear- 
ance of the QPO is not a consequence of the transition, 
but possibly its cause), while the QPO frequency decreases 
monotonously; the color radius first decreases before it 
starts growing, when the count rate decreases, marking 
the transition to the low state: this causes the reverse cor- 
relation seen in Fig. |8| On the other hand, returning to 
the data points in the high state, before the onset of the 
QPO one sees that in its variations the color radius consis- 
tently shows a fixed lower bound, which probably marks 
the Last Stable Orbit. 

The data points showing an inverse QPO frequency-color 
radius correlation lie at lower color radius, and at these 
times the disk temperature decreases, as predicted by 
MFR. It is thus very likely that, in the spectral fits for 
these points, the correlation between the color radius and 
the disk inner radius has become anomalous. 



3. Discussion 

The possibility, which will be explored in Paper II, of a 
theoretical explanation has led us to reconsider the con- 
trasting correlations, observed by SMR in GRO J 1655-40 
and XTE J1550-564, between the disk color radius and the 
QPO frequency First guided by the arguments of MFR, 
we have rejected data points for GRO J1655-40 where 
the measure of the disk color radius was dubious. The 
major result of this work is thus that, though we have 
used the most stringent criterion to select "good" obser- 
vational data points, the reverse correlation remains: the 
QPO frequency decreases with decreasing radius in the 
case of GRO J1655-40. 

In the case of XTE J1550-564 the data, for the observa- 
tions where a QPO is seen, always correspond to a low 
value of the ratio disk black body/total flux, and may not 
be reliable according to the criterion of MFR. We have 
thus turned to GRS 1915+105, because previous work by 
SM97 and Muno et al. (1999) had shown that, during the 
frequently observed 30 min. cycles of that source, the QPO 
frequency and disk color radius varied in a correlated man- 
ner over a large interval. We have done a complete analy- 
sis of two of these cycles, sampling data over the minimal 
time interval compatible with proper statistics in order 
to follow as precisely as possible the rapid variations of 
the source. This does not show a reverse correlation as in 
GRO J 1655-40. Data points obtained at the onset of the 
QPO would seem to indicate it, but a detailed analysis 
shows that they are probably incorrect (the spectral fit 
returning an anomalous value of the disk color radius), so 
that we have to rule them out. 

However the measure of the radius is not an absolute 
one. It is obtained from a multicolor black-body + power 



Rodriguez, Varniere, Tagger & Durouchoux: AEI & QPO in BH Binaries 




0.0100 
0.0010 



0.1 1 10 100 0.1 1 10 100 0.1 1 10 100 0.1 1 10 100 0.1 1 10 100 0.1 1 10 10 



Fig. 6. Power Spectra from the second cycle; top left corresponds to the first time interval analyzed in this cycle, 
and bottom right to the last one. Arrows indicate the position of the QPO found from the fit in each interval, using 
the dynamical power spectrum (figure ^|) as a guide to reduce ambiguities. Even with this help, the limited statistics 
makes it difficult or impossible in some cases to correctly identify the QPO. 
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Fig. 7. Plot of the timing evolution of the source flux during the first (left), and second (right) cycle. Subpanel is a 
temporal zoom on the interesting timing region, with, under the zoomed lightcurve, the fit values of the QPO frequency 
and the disk radius. Error bars are one a confidence level. 
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Fig. 8. Plot of vqpo vs r co i for GRS1915+105 for the first 
(squares) and second (circles) cycles; error bars shown on 
figure fj] arc not displayed here, for the sake of clarity. 

law tail model (plus an iron line when needed), involving 
a number of hypotheses (such as axisymmctry of the disk 
emission, and ~ r~ 3 / 4 radial temperature profiles), and 
should be submitted to a number of corrections (such as 
electron diffusion or spectral hardening factor) which we 
have not tried to include here. Furthermore, MFR have 
recently shown that effects due to the finite disk thickness 
could result in a systematic underestimation of the radius, 
when it becomes small. It is remarkable here that, for the 
GRS 1915+105 data we have rejected, the spectral fit re- 
turns a color temperature which is lower than usual, and 
in particular lower than at times just preceding the onset 
of the QPO, although the radius decreases: this is compat- 
ible with the predictions of MFR, whose mechanism relies 
in part on the dissipation of part of the accretion energy 
in the corona rather than the disk (so that the disk cools 
down) . 

On the other hand, the data we have rejected for GRO 
J1655-40 show a very different behavior, also frequently 
observed in the low and steady state of GRS 1915+105 
(Muno et ai, 1999): there, the color radius takes ex- 
tremely small values (sometimes very few kilometers) 
while the color temperature becomes high. This is not 
the behavior predicted by MFR, and it is quite unlikely 
that corrections of this type could explain so strongly 
anomalous values of the radius: the lowest radius value 
found by Sobczak et al. (2000) is about 4 times lower 
than the ones we retain. 

However, we note that, if the emissivity is close to 
unity, the spectral fits actually measure (essentially by 
the intensity of the black body emission) the size of the 
emitting region: a simplistic view of the multicolor black- 
body spectrum used in the spectral fits is that the peak 
of the emission corresponds to the highest temperature in 
the disk, and the total intensity to the area of the emit- 
ting region at that temperature, all this weighted by the 
temperature distribution in the disk. In GRO J1655-40 we 
find that the color radius can strongly vary, while the to- 
tal luminosity remains approximately constant. The most 
obvious explanation would then be that the total accre- 
tion energy dissipated in the disk does not vary much, 
but that the area of the region where it is dissipated does. 
This leads us to suggest that the anomalously low mea- 
sured radius, together with a high temperature, could be 
interpreted by the presence in the disk of a hot point or a 
spiral shock, where a substantial fraction of the accretion 
energy would be dissipated. A spiral shock could result 
from the non-linear evolution of the AEI, just as the gas 
forms shocks in the arms of spiral galaxies. 
The discussion about the exact ratio r co i/ri n , related to 



the hardening factor by r in = r co/ x / 2 , with / = T col /T e ff 
, being the hardening factor (Ebisawa et ai, 1994) is a very 
complex one. It appears to be dependent on many phys- 
ical parameters, as shown e.g. by Shimura & Takahara 
(1995), for a Schwarschild Black Hole. Although both GRS 
1915+105 and GRO J1655-40 are thought to be a Kerr 
Black Hole, we adopted the best result found in their work, 
/ = 1.9, and tentatively estimated the position of the last 
stable orbit, using Rlso ~ § x 0.8i? co ;/ 2 (Ebisawa et ai, 
1994). We chose for each source the lowest "good" value of 
R co i as the one approaching Rlso the most. In the case of 
GRO J1655-40, this leads to an estimation of Rlso ~ 36 
km, inconsistent with a ~ 7M Q Schwarschild Black Hole, 
where Rlso would be of the order of 63 km. This large 
difference can be due to two effects. First, the value we 
adopt is not the exact measure of the last marginally sta- 
ble orbit, but is really a rough estimation of it (we at least 
expect it to be closer to the compact object). The second 
one, which may be the major effect, is that this source, in 
particular, is expected to be a maximally rotating Black 
Hole (Strohmaycr, 2001), where Rlso = 17 km, which is, 
then, much more consistent (Paper 2). The case of GRS 
1915+105 can now be explored, since Greiner (2001) has 
estimated the mass of the Black Hole to be close to 14M©. 
The result found {Rlso ~ 59 km) , as for GRO J1655-40, 
is inconsistent with a Schwarschild Black Hole of 14M©, 
but would perfectly fit the case of a maximally rotating 
Kerr Black Hole of such mass (see Paper 2). We also 
note as an interesting possibility to discuss the mass of 
the Black Hole the work of di Matteo and Psaltis (1999), 
which should be enriched by recent measurements of this 
mass. 

But we have to exercise here, extreme caution, since : 

— all the calculations in Shimura & Takahara (1995), in 
particular the determination of the best value of the 
hardening factor, are made for a Schwarschild black 
hole; 

— all the calculations, presented in Ebisawa et al. (1994), 
and Shimura & Takahara (1995), are done for a 
Shakura & Sunyaev (1973) a-disk, where energy and 
angular momentum are locally dissipated in the disk 
under the effect of viscous stresses; whereas in the 
framework of the AEI, only a small part of the en- 
ergy warms the disk. Most of it is transported toward 
the corotation radius (Tagger & Pellat, 1999; Paper 2, 
Varniere et al., in Prep.); 

— the Kerr metric may also affect the relation of Ebisawa 
et al. (1994), between Rlso, the radius of the Last 
Stable Orbit, and the measured R co i', 

— although in the case of GRO J 1655-40 the black hole is 
thought to be rapidly rotating, this is still unsolved for 
GRS 1915+105; its temporal behavior does not allow 
a firm conclusion (Strohmaycr, 2001) (one has to note, 
nevertheless, that previous studies give preference to a 
rapidly rotating black hole (Zhang et al. , 1997). 
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The relevance of these calculations is thus a matter 
of debate, and we cannot draw conclusion on the value 
/ = 1.9 for the hardening factor, or on the rotation of 
the black holes. It is however remarkable that the inde- 
pendent method presented in Paper 2 gives a good range 
of masses, for both sources, in the cases where both are 
spinning rapidly. 

Our conclusions must be considered as tentative, given the 
uncertainties on the measure of the disk color radius and 
its relation with its inner radius, and given the other ef- 
fects associated with the disk parameters, such as the am- 
plitude and the radial distribution of the magnetic field. 
On the other hand, a coupled investigation of the low fre- 
quency QPO and of the disk spectral properties seems to 
be a most promising way to probe the physics of the inner 
accretion disk in black hole binaries. 
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